Zero-voltage conductance peak from weak antilocalization in a Majorana nanowire 



D. I. Pikulin, 1 J. P. Dahlhaus, 1 M. Wimmer, 1 H. Schomerus, 2 and C. W. J. Beenakker 1 

1 Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 
2 Department of Physics, Lancaster University, Lancaster, LAI 4YB, United Kingdom 

(Dated: August 2012) 

We show that weak antilocalization by disorder competes with resonant Andreev reflection from a 
Majorana zero-mode to produce a zero- voltage conductance peak of order e 2 /h in a superconducting 
nanowire. The phase conjugation needed for quantum interference to survive a disorder average is 
provided by particle-hole symmetry — in the absence of time-reversal symmetry and without requir- 
ing a topologically nontrivial phase. We identify methods to distinguish the Majorana resonance 
from the weak antilocalization effect. 
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I. INTRODUCTION 

Weak localization (or antilocalization) is the system- 
atic constructive (or destructive) interference of phase 
conjugate series of scattering events. In disordered met- 
als it is time-reversal symmetry that provides for phase 
conjugation of backscattered electrons and protects their 
interference from averaging out to zero [TJ [5]. A mag- 
netic field breaks time-reversal symmetry, changing the 
disorder-averaged conductance by an amount SG of or- 
der e 2 /h. The sign of SG distinguishes weak localization 
(SG < 0, conductance dip) from weak antilocalization 
(SG > 0, conductance peak). 

Andreev reflection at a superconductor provides an 
alternative mechanism for phase conjugation due to 
particle-hole symmetry. No time-reversal symmetry is 
needed, so weak (anti)localization can coexist with a 
magnetic field and is only destroyed by a bias voltage 
[31 13] ■ The resulting zero-bias anomaly in the conduc- 
tance of a normal-metal-superconductor (NS) junction 
is obscured in zero magnetic field by the much larger ef- 
fects of induced superconductivity, which scale with the 
number of transverse modes N in the junction. These 
order Ne 2 /h effects are suppressed by a magnetic field, 
only the order e 2 /h effect from weak (anti)localization 
remains [5]. 

In a superconducting nanowire there is an altogether 
different origin of zero-bias anomalies in a magnetic field, 
namely the midgap state that appears at the NS inter- 
face following a topological phase transition [SHE] ■ Reso- 
nant Andreev reflection from the zero- mode gives a 2e 2 /h 
conductance peak at zero voltage [5]- The first reports 
[T0HT2] of this signature of a Majorana fermion are gen- 
erating much excitement |13| . There is an urgent need 
to understand the effects of disorder, in order to deter- 
mine whether it may produce low-lying resonances that 
obscure the Majorana resonance [T4UT8] . 

These recent developments have motivated us to in- 
vestigate the interplay of Majorana zero-modes and 
weak (anti)localization. Earlier studies of weak 
(anti) localization at an NS junction [3] l4l HMSTj did not 
consider the possibility of a topologically nontrivial phase 



with Majorana fermions. Calculations of the local den- 
sity of states near a zero-mode [22H24] address the same 
physics of midgap quantum interference that we do, but 
cannot determine the conductance. 

This paper consists of two parts: We first give a simple 
model of a disordered NS interface that allows us to ob- 
tain analytical results for SG with and without Majorana 
zero-modes. We then turn to a numerical simulation of a 
Majorana nanowire and compare the conductance peak 
due to weak antilocalization (in the topologically triv- 
ial phase) with that from a Majorana zero-mode (in the 
nontrivial phase). The two effects can appear strikingly 
similar, but we will discuss several ways in which they 
may be distinguished. 

Before we present our findings, we wish to emphasise 
that it is not the purpose of this work to diminish the sig- 
nificance of experiments reporting the discovery of Ma- 
jorana fermions in superconductors. On the contrary, we 
feel that existing [UH - tTS] and forthcoming experiments 
will gain in significance if possible alternative mecha- 
nisms for zero-voltage conductance peaks in a magnetic 
field arc identified and understood, so that they can be 
ruled out. Weak antilocalization was so far overlooked as 
one such mechanism. 



II. ANALYTICAL THEORY 

For the analytical calculation we consider a super- 
conducting wire that supports Q topologically protected 
zero-modes at the interface with a normal metal (see 
Fig. [IJ . We assume broken time-reversal symmetry (by 
a magnetic field) and broken spin-rotation symmetry (by 
spin-orbit coupling), so generically |25) this is symme- 
try class D with Q £ {0,1}. If the wire is sufficiently 
narrow (relative to the spin-orbit coupling length), an 
approximate chiral symmetry [551 HZ] promotes the su- 
perconductor to class BDI with Q e {0, 1, 2, . . . N}. The 
integer N is the number of propagating electronic modes 
through the wire in the normal state, counting both spin 
and orbital degrees of freedom. 
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The matrices U, V, U', V are NxN unitary matrices and 
W T = diag (Ti , T2 , . . . Tjv) is a diagonal matrix of transmis- 

sion eigenvalues of the normal region. In class BDI chiral 
symmetry relates V — U T , V = V T . 



FIG. 1: A bias voltage V applied to the normal metal (N) 
drives a current / into the grounded superconductor (S). Elec- 
trons and holes (e, h) are scattered by disorder or a tunnel 
barrier in N and converted into each other by Andreev re- 
flection at the NS interface, as described by the scattering 
matrices 5n and rA- Particle- hole symmetry ensures that the 
phase shifts accumulated by e and h along a closed trajectory 
cancel, irrespective of whether time-reversal symmetry is bro- 
ken or not. Such phase conjugate series of scattering events 
permit weak (anti)localization to persist in a magnetic field. 



A. Scattering matrix 

We construct the scattering matrix of the NS junction 
at the Fermi level by assuming a spatial separation of nor- 
mal scattering in N and Andreev reflection in S. Within 
the excitation gap there is no transmission through the 
superconductor. The matrix of Andreev reflection 
amplitudes from the superconductor is then a 2N x 2N 
unitary matrix. Mode mixing at the NS interface can be 
incorporated in the scattering matrix Sn of the normal 
region, so we need not include it in rx- It has the block 

form nana 



M , 

'■■>=i a- r =e( c T 

m=l V 





cos a r , 





i sin a Tl 



—1 sin a r , 




(1) 



The symbols l n ,0n denote, respectively, an n x n unit 
matrix or null matrix. We have defined f = if the 
difference N — Q is even and £ = 1 if N ~ Q is odd, so 
that N — Q — ( = 2M is an even integer. The Andreev 
reflection eigenvalues p m = sin 2 a m that are not pinned 
at or 1 are twofold degenerate |50) . 

The normal region has scattering matrix 
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The electron and hole blocks (with NxN reflection and 
transmission matrices r,r',t,t') are each others complex 
conjugate at the Fermi level. The off-diagonal blocks of 
5n vanish, because the normal metal cannot mix elec- 
trons and holes. The matrix so is unitary, sqSq = 1, 
without further restrictions in class D. In class BDI chi- 
ral symmetry requires that sq = sj is also a symmetric 
matrix. 

To separate the mixing of modes from backscattering, 
we make use of the polar decomposition 
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B. Conductance 

We combine Sn and ta to obtain the matrix r^ e of 
Andreev reflection amplitudes (from electron e to hole h) 
of the entire system. This calculation is much simplified 
in the case £ = 0, p m = 1 (m = 1,2, ... M) that all 
modes at the NS interface are Andreev reflected with 
unit probability. For this case T = 0, N — Q = 2M, we 
obtain 



The 



t'*A*(l-rAr*A*)-H, A — cr® M © 1c 



(4) 



The notation c® M signifies the 2M x 2M matrix con- 
structed as the direct sum of M Pauli matrices. 

The Andreev reflection matrix determines the conduc- 
tance 



G = G Trr he r{ e , G Q = 2e 2 /h. 



(5) 



Substitution of the polar decomposition ([3| gives the 
compact expression 



G/G = Tr TMTM\ 

M = {i-n*VT^TfiVT zr T~r 1 n*, n = v'Av* 



(6) 



This is the conductance at the Fermi level, in the limit 
of zero bias voltage. Away from the Fermi level particle- 
hole symmetry is broken, so the electron and hole blocks 
in Sn are distinct unitary matrices s e and s^. If the 
voltage remains small compared to the excitation gap, 
we can keep the same r\. The finite- voltage differential 
conductance G = dl/dV is then given by 



G/G = Tr T h MT e M\ 

m = (i- niVi-Te^Vi-Th)- 1 ^, 
^ = M, n h = v,[Av:. 



(7) 



Chiral symmetry remains operative away from the Fermi 

flh = fit in class BDI. 



level, hence V' e = Vj, V } [ 



C. Random matrix average 

Isotropic mixing of the modes by scattering in the nor- 
mal region means that the unitary matrices in the po- 
lar decomposition ^ are uniformly distributed in the 
unitary group U(N). We can calculate the average con- 
ductance for a given set of transmission eigenvalues by 
integration over U(N) with the uniform (Haar) measure. 
A full average would then still require an average over 
the TVs, but if these are dominated by a tunnel barrier 
they will fluctuate little and the partial average over the 
unitary matrices is already informative. 
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FIG. 2: Amplitude SG of the average zero-voltage conduc- 
tance peak as a function of (mode-independent) transmission 
probability T, in symmetry class D (thick curves) and BDI 
(thin curves) for different number of modes N. The super- 
conductor is topologically trivial when N is even (Q — 0, 
dashed curves) and nontrivial when N is odd (Q = 1, solid 
curves). The dash-dotted curve is the Q-independent large- N 
limit p3| . 

The calculation is easiest if all T n 's have the same value 
< T < 1. The average zero- voltage conductance (G) is 
then given by the integral 

(G)=T 2 G d<t> P {<p)\l-(\-Ty\\ (8) 
Jo 

with p(4>) = (Y]„ S(<p — (pn)) the density on the unit circle 
of the eigenvalues e^" of ftft* . The corresponding finite- 
voltage expression has a uniform p = N/2tt, leading to 

(G) = NGqT/ (2 — T), (9) 

irrespective of the symmetry class and independent of the 
topological quantum number Q. 

The zero-voltage average (|8| does depend on Q and is 
different for class D and BDI. The calculations are given 
in the Appendix. Explicit expressions in class D are 

{2T for N = 2, Q = 0, 

1 + 2T 2 for N = 3, Q = 1, 

2T{2-T + T 2 ) for N = A, Q = 0, 

1 + 2T 2 (3-2T + T 2 ) for A = 5, Q = l. 

(10) 

The Q-dependence appears to second order in the reflec- 
tion probability R = 1 — T, while the general first-order 
result 

(G/G ) D = N(1-2R) + 2R + 0(R 2 ) (11) 

is Q-independent. The corresponding expressions in class 
BDI are more lengthy, and we only record the small- R 
result 

O 2 4- N 

(G/G )bdi = N(l - 2R) + 2R^— \— + 0(R 2 ), (12) 



to show that it is Q-dependent already to first order in 
R. These are all finite-AT results. In the large-A^ limit 
the Q-dependence is lost, 

NT 9(^ — T) 
(G/G ) = — + + OiN- 1 ), (13) 

irrespective of the symmetry class. 

As illustrated in Fig. [2] for this case that all T n 's have 
the same value T the difference SG = (G) — (G) is pos- 
itive, corresponding to weak aniilocalization and a con- 
ductance peak. The sign of the effect may change if the 
T„'s are very different, in particular in class BDI — which 
has SG < in a quantum dot geometry (circular ensem- 
ble) [H] . This is a special feature of quantum interference 
in a magnetic field, that the distinction between weak lo- 
calization and antilocalization is not uniquely determined 
by the symmetry class [3TJ 132 E2] • 

III. SIMULATION OF A MICROSCOPIC 
MODEL 

A. Model Hamiltonian 

We now turn to a microscopic model of a Majorana 
nanowire [7] : a conducting channel parallel to the x- 
axis on a substrate in the x-y plane (width W , Fermi 
energy E-p), in a magnetic field B (orientation h, Zee- 
man energy Ez — ^g e sPBB), with Rashba spin-orbit 
coupling (characteristic energy E so — m c sa 2 a /h 2 , length 
Z so = h 2 /m e ffa so ), and induced s-wave superconductivity 
(excitation gap Ao). The Hamiltonian is 

V _(H - E F Aa y \ 
{ AV, E F -HZ)> 

H = ^ — + U(x,y) + ^-((TxPy - (TyPx) + E z h ■ a. 

(14) 

The electrostatic potential U = J7 ga t G + SU contains 
the gate potential t/ ga t that creates the tunnel barrier 
and the impurity potential SU that varies randomly from 
site to site on a square lattice (lattice constant a), dis- 
tributed uniformly in the interval (— [/disorder, ^disorder)- 
The disordered region is — Ln < a; < Lg, an NS interface 
is constructed by increasing the pair potential A from 
to Ao at x = 0, and a rectangular barrier of height 

[/"barrier, thickness <5L b arrier, IS placed at X = -^barrier- 

The conductance of the normal region (x < 0) contains 
a contribution Gdisordcr from disorder and Gbamcr from 
the barrier. 

The orientation of the magnetic field plays an impor- 
tant role [3 [7]: It lies in the x-y plane to eliminate orbital 
effects on the superconductor and we will only include its 
effect on the electron spin (through the Zeeman energy) . 
A topologically nontrivial phase needs a nonzero excita- 
tion gap for Ez > Ao, which requires a parallel magnetic 
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FIG. 3: Disorder-averaged differential conductance as a function of bias voltage, for a nanowire modeled by the Hamiltonian 
(14 1. The two panels a) and b) correspond to the two geometries shown to scale above each plot. (The solid vertical line 
indicates the position of the tunnel barrier, relative to the NS interface; disordered regions are dotted.) Each panel shows data 
for zero magnetic field (black), and for two nonzero magnetic field values (blue and red). The solid curves are for parallel field 
Bu and the dashed curves for perpendicular field B±. The system is topologically trivial (Q — 0) in all cases except for the red 
solid curves (Q = 1). 



field Bn (h — x). We will consider that case in the next 
subsection, and then discuss the case of a perpendicular 
magnetic field B± [n = if) in Sec. Ill C 



B. Average vs. sample-specific conductance 

To avoid the complications from chiral symmetry we 
first focus on a relatively wide junction, W = 3l S o, when 
symmetry class D (rather than BDI) applies 29 . (We 
turn to class BDI in the next subsection.) The normal 
region has N = 8 propagating modes (including spin) 
in zero magnetic field, for Ep = 12E so . The topologi- 
cal quantum number Q was determined both from the 
determinant of the reflection matrix |331 134] , and inde- 
pendently by counting the gap closings and reopenings 
upon increasing the magnetic field. A transition from 
Q = to Q = 1 is realized by increasing B\\ at fixed 
A = 8 E so . 

Results are shown in Fig. [3] (solid curves) for two ge- 
ometries, one with the tunnel barrier far from the NS 
interface (^barrier = 23 Z so , SL halIlcI = 4a = 0.4 Z so , 

^disorder 12.5 E so ^ U barrier 1^ E SQ7 ^disorder 

2.7e 2 /h 1 Gbarricr = l.Se 2 /h) and another with the bar- 



— 3 Z S o; barrier & 
50 E SQl Gdisordcr 



rier close to the interface (xh& 

0-1 ^soj ^disorder — 20 E SOl Uh, 

2.9e 2 /h, G barricr = 3.2e 2 //i). 

The disorder-averaged conductance shows a zero- 
voltage peak in a magnetic field, regardless of whether 
the nanowire is topologically trivial (Q = 0) or nontriv- 
ial (Q = 1). The peak disappears in zero magnetic field 
and instead a conductance minimum develops, indicative 
of an induced superconducting minigap in the normal re- 
gion. The two geometries in panels [3^, andJSJa show com- 



parable results, the main difference being a broadening 
of the zero-bias peak when the tunnel barrier is brought 
closer to the NS interface — as expected from the in- 
crease in Thouless energy |35) . The shallow maximum 
which develops around zero voltage in the B — curve 
of panel [3]d is a precursor of the reflectionless tunneling 
peak, which appears in full strength when the barrier is 
placed at the NS interface [5]. 

This all applies to the average conductance in an en- 
semble of disordered nanowires. Individual members of 
the ensemble show mesoscopic, sample-specific conduc- 
tance fluctuations, in addition to the systematic weak 
antilocalization effect. For some disorder realizations the 
zero-voltage conductance peak remains clearly visible, 
see Fig. [4] The peak sticks to zero bias voltage over 
a relatively wide magnetic field range, even though the 
superconductor is topologically trivial (Q = 0). The ap- 
pearance and disappearance of the peak is not associated 
with the closing and reopening of an excitation gap, so 
it cannot produce Majorana fermions [36] . 



C. Parallel vs. perpendicular magnetic field 

So far we considered a class D nanowire with magnetic 
field Bn parallel to the wire axis. In a perpendicular mag- 
netic field B± (perpendicular to the wire in the plane 
of the substrate) the symmetry class remains D (bro- 
ken time-reversal and spin-rotation symmetry) , although 
the topologically nontrivial phase disappears [6l [7] . We 
therefore expect the class D zero-bias peak to persist in a 
perpendicular field as a result of the weak antilocalization 
effect. 

This expectation is borne out by the computer sim- 



5 





0.04 


u 


0.02 









u 








-0.02 




-0.04 









b) 



0.4 0.8 1.2 

Ez [E so ] 



1 2 3 
Ez [E so ] 



0.4 



0.2 





-0.2 
-0.4 



FIG. 4: Numerical simulation of a nanowire for a single disorder realization (no averaging). The color scale gives the differential 
conductance as a function of bias voltage (vertical axis) and parallel magnetic field (horizontal axis). The parameters in panels 
a,b correspond to those in Fig. [3^,b. The magnetic field range in both panels is in the topologically trivial phase (Q = 0), but 
still exhibits a conductance peak pinned to zero voltage (green circle). 
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FIG. 5: Same as Fig. |3j but now for a narrower wire in sym- 
metry class BDI (rather than D). The system is topologically 
trivial, without Majorana zero-modes. The weak antilocal- 
ization peak vanishes if the magnetic field is rotated from Bu 
to B±. (The parameter values are listed in Ref. [37] . 'i 



ulations, see the dashed curves in Fig. [3] A zero-bias 
peak exists for both B± and B» . If the nanowire is topo- 
logically trivial, there is not much difference in the peak 
height for the two magnetic field directions (compare blue 
solid and dashed curves). In contrast, if the nanowire is 
topologically nontrivial for parallel field then the peak 
is much reduced in perpendicular field (red solid versus 
dashed curves). The disappearance of the Majorana zero- 
mode and the collapse of the zero-bias peak may be ac- 



companied by the appearance of propagating modes in 
the superconducting part of the nanowire. This explains 
the increased background conductance in the red dashed 
curve of Fig. [3^,. 

The effect of a magnetic field rotation is entirely differ- 
ent when W < l so and the symmetry class is BDI rather 
than D [H] . The term a x p y in the Hamiltonian ( 14 ) 
can then be neglected, so that W commutes with a y in 
a perpendicular magnetic field (n = y). The two spin 
components along ±y decouple and for each spin com- 
ponent separately the particle-hole symmetry is broken. 
We therefore expect both the Majorana resonance and 
the weak antilocalization peak to disappear in a perpen- 
dicular magnetic field for sufficiently narrow wires. 

This is demonstrated by the computer simulations 
shown in Fig. [5j for the average conductance in a topo- 
logically trivial wire of width W = 0.3 l so . The main dif- 
ference with the data in Fig. [3] is that the symmetry class 
is now BDI rather than D, because of the narrower wire. 
This change of symmetry class does not significantly af- 
fect the weak antilocalization peak in a parallel magnetic 
field. But if the magnetic field is rotated to a perpendic- 
ular direction, the peak disappears — as expected for a 
class BDI nanowire. 



IV. DISCUSSION 

In conclusion, we have shown that random quantum in- 
terference by disorder in a superconducting nanowire can 
systematically produce a zero-voltage conductance peak 
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FIG. 6: Differential conductance for a single disorder realization of a nanowire (N = 2 spin-resolved modes, parameter values 
are listed in Ref. [39]). The left panel shows the appearance of a zero-voltage peak in a range of magnetic field values, for B 
parallel to the wire. The right panel shows the dependence on the orientation of the magnetic field, for a fixed field strength 
(Ez = 10 -E so ). The zero- voltage peak vanishes if B is perpendicular to the wire. This is the same phenomenology as for a 
Majorana resonance, but here it happens in the topologically trivial phase. 



in the absence of time-reversal symmetry. This weak an- 
tilocalization effect relics on the same particlc-hole sym- 
metry that protects the Majorana zero-mode, but it ex- 
ists in both the topologically trivial and nontrivial phase 
of the superconductor. A conclusive demonstration of 
Majorana fermions will need to rule out this alternative 
mechanism for a conductance peak. 

There are several strategies one might follow for this 
purpose: 

• Increasing the tunnel barrier with a gate voltage 
suppresses the weak antilocalization effect, but not 
the Majorana resonance. The resonance does be- 
come narrower, so at finite temperatures thermal 
smearing will still lead to a suppression with in- 
creasing barrier height and this might not be the 
most effective strategy to distinguish the two ef- 
fects. 

• The disappearance of the conductance peak when 
the magnetic field is rotated (in the plane of the 
substrate) towards a direction perpendicular to the 
wire, the technique used in Refs. [TOIIH], can iden- 
tify the Majorana zero-mode — but only if the ratio 
W/l so is sufficiently large that the wire is in class D 
rather than BDI. In class BDI the Zeeman energy 
in the rotated field commutes with the Rashba en- 
ergy, precluding the weak antilocalization effect as 
well as the Majorana resonance. Both Refs. [TOj, [12] 
have W < l so and are believed to be in class BDI 
[161 126] . so this complication seems quite relevant. 



• Measuring the conductance through a single-mode 
point contact is a very effective strategy: for N = 1 
the zero-temperature conductance G — Q x 2e 2 /h 
directly measures the topological quantum number 
even without any tunnel barrier [38], and this sig- 
nature of a Majorana zero-mode is quite robust 
against finite temperatures. The single mode in 
the point contact should be spin resolved for this 
to work: If instead the point contact transmits both 
spins in one orbital mode (N — 2), then the ambi- 
guity between weak antilocalization and the Majo- 
rana resonance remains (see Fig. KM . 
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Appendix A: Random-matrix theory 

To evaluate the average conductance Q we seek the 
density of the eigenvalues x n = e l< ™ of the product 
X = of the unitary matrix f2 and its complex con- 
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jugate. We denote \i n = cos0„ e [—1, 1] and determine 
the joint probability distribution P({/i„}) using methods 
from random-matrix theory [40] . 

In symmetry class D, we have f2 = V'AV* with V and 
V independently and uniformly distributed according to 
the Haar measure dU of the unitary group U(N). Be- 
cause d(UU') = dU for a fixed unitary matrix U', the 
matrix = U is itself uniformly distributed in U(N). 

In class BDI, we have V = V T and we may write 
il = UXW with U uniformly in U(N). The diagonal 
matrix A = diag(Ai, A2, • • • A at) contains the eigenvalues 
A n = ±f of A. The number q — \Q\ of Majorana zero- 
modes is encoded in the topological invariant Q — Tr A = 
J2 n ■ (For full generality we allow Q to also take on 
negative values, but the final result will only depend on 
the absolute value q.) 



1. Brownian motion of unitary matrices 

We employ Dyson's Brownian motion approach [41], 
which sets up a stochastic process for the unitary ma- 
trix U whose stationary distribution coincides with the 
Haar measure on U{N). In each infinitesimal step of the 
process, U —¥ U exp(iH), where H is a Hermitian matrix 
from the Gaussian unitary ensemble, with identically nor- 
mal distributed complex numbers Hi m = B* ml (I < to), 



conjugation is basis dependent; if B = A* in one ba- 
sis then this relation is only preserved under orthogonal 
transformations, but not under general unitary trans- 
formations. Thus, we work in a fixed basis |r) (at 
most permitting orthogonal basis changes), and define 
for any \ip) — J2 r ^r\f) a complex-conjugated vector 
W) = Er#k>. As usuaL M = ErV£<r|; thus 
(r\=E r Mr\. 

The matrices X and U are unitary and obey Det X = 
|DetC/| 2 = 1. Moreover, the matrix X* has the same 
eigenvalues x%, 22, . . . scy as the matrix X. For even N, it 
follows that all eigenvalues appear in complex-conjugated 
pairs; every eigenvalue x k has a partner xj, — x* k = x% . 
For odd N, in addition to such pairs there is a single 
unpaired eigenvalue, denoted as %, which (because of 
the constraint on the determinant) is pinned at xn = 1. 
The paired eigenvectors are related according to 



\k) = Z k U\k* 



(A2) 



Here we have to set such that £| = Afe; this guarantees 
that the relation between both eigenvectors in a pair is 
reciprocal, \k) — \k). Observing that the eigenvectors 
form an orthogonal basis, we find the matrix elements 



(k\u\n = tkS k i= ((k*\u*\i)y = (i\u T \k* 



(A3) 



0, HkiH v 



i 



SkmSinT; the 



With help of these matrix elements we can now eval- 
uate the drift and diffusion coefficients. In second-order 
perturbation theory, 



is implied to generate infinitesimal incre- 



5xi = (l\SX\l) 



X[ - x k 

where the prime restricts the sum to k 7^ / while 



(A4) 



Hi, 

limit r 
ments. 

The corresponding increments <5/i n can be calcu- 
lated in perturbation theory. The drift coefficients 
c; = lim r _>.o T _1 <5/i/ and the diffusion coefficients c/ m = 
lim T ^ T^SmS/im follow by averaging over the random SX = iUHU* - iXH* + UHU*H* - \UH 2 U* - \XE 



y2>{l\5X\k){k\5X\l) 
k 



variables in H. As we will see, the symmetries in the 
classes D and BDI are restrictive enough so that these 
coefficients can be expressed in terms of the quantities 
/i„ alone, without requiring data from the eigenvectors 
of X. Thus, the stochastic process for these quantities 
closes. 

Introducing a fictitious time t, the evolution of the joint 
probability distribution is governed by a Fokker-Planck 
equation, 



dP 
~dt 



E 



_d_ 
dm 



ci + 



d d 



2 ^ d[ii d[i, 



P({»n},t) 



(Al) 

The stationary solution P({// n }), for which the right- 
hand-side of the Fokker-Planck equation vanishes, is the 
required eigenvalue distribution. 



2. Symmetry class D 

In class D we have X = UU* with U uniformly dis- 
tributed in U (N) . Notice that the operation of complex 



(A5) 

is the increment of X to leading order in r. The Gaussian 
averages are now carried out according to the rules 

(k\AHB\l)(m\CHD\n) = r(k\AD\n}(m\CB\l), (A6) 

(k\AHB\l)(m\CH*D\n) = T(k\AC T \m*)(n*\D T B\l) . 

(A7) 

In particular, IP = Nt, UHU*H* = tUW = t, and 

(l\UHU* - XH*\k}(k\UHU* - XH*\l) 

= 2r(l\X\l)(k\X\k) - T(l\U T \k*)(l*\U*\k) 

-T(k\U T \l*)(k*\U*\l) 
= 2rxix k - rS lk -(xi + xj), (A8) 

where we invoked Eq. ( |A3[ ). We thus obtain 
Ix-^r-Nrx^rY!- 1 ^' 61 ^^^ 

k 

Analogously, we find 



xi - x k 



(A9) 



6xi5x m = (l\SX\l) (m\6X\m) = —2r5i m x l +2rSim. 



(A10) 



Note that these expressions only depend on the eigen- 
values. We remark that for the pinned unpaired eigen- 
value x n = 1, occurring if N odd, these relations deliver 

$%n — {3xn) 2 = 0. 

We now pass over to the quantities /Ltj = (xi + x{)/2, 
and restrict the index I such that it enumerates the pairs 
of eigenvalues. For even N we then find 



1 



Mi - Mfc 



while for odd N we have 



-3r W - 2t(m? - 1)^' 



Mi - Mfc 



(All) 



(A12) 



where the double-prime excludes the pinned eigenvalue. 
Furthermore, 



5fj,iSfx m = 2r(l - fifydit, 



(A13) 



The stationarity condition of the associated Fokker- 
Planck equation (Al) can be expressed as 
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1 d 



dfn 2 dm 

For even N = 2M, this is solved by 

M 



(A14) 



P(/ii,/i 2 , ■ • - Mm) oc 



M M 

II. II » 



k=l \A ~ Mfc I. ,,,-] 

(A15a) 

up to a normalization constant. Each of the Mn' s (n = 
1,2,... A/) is twofold degenerate. For odd N = 2M + 1 
one eigenvalue is pinned at +1, and the remaining ones 
are twofold degenerate with distribution 



M 



M 



H — UrnY 



P(lH,fJ>2, •••Mm) oc JJyl-Mfc II (M/ 

k=l ;<rn=l 

(A15b) 

This concludes our derivation of the eigenvalue distri- 
bution of UU* with U uniform in U(N). We have not 
found the result (A15) in the literature, but there is a 



curious correspondence with the known |40[ 142] eigen- 
value distribution of orthogonal matrices (uniformly dis- 
tributed according to the Haar measure). An (N + 1) x 
(N + 1) orthogonal matrix O with determinant —1 has 
one eigenvalue pinned at —1. If we exclude that eigen- 
value, the remaining N eigenvalues of O have same prob- 
ability distribution as the N eigenvalues of UU* . 



The Brownian motion is induced by 0(1 + A + A 2 /2), 
where (in the fixed basis) A = —A T is a real antisym- 
metric matrix, with Af m — r. Due to the condition on 
the determinant, there is always one eigenvalue pinned at 
xn+i = — 1, while an additional eigenvalue is pinned at 
xn = 1 if N is odd. All other eigenvalues appear in pairs 
xi, xj, with \l) = \l*) (no additional factors are required). 

We calcaluate the increments and average: 



8 Xl = \{l\OA 2 \l)+Y. 



Sxi = — | tNxi + t 



(l\OA\k)(k\OA\l) 
xi - x k 

xix k (S kT -l) 



k^l 



xi - x k 



(A16) 



SxiSx k = (l\OA\l)(k\OA\k) 



=> Sx t 5x k = Txix k {8^ - 8i k ) = t{5 1 j, - x t Si k ). (A17) 

(Note that (Z|A|i) does not vanish if \l) is complex, as is 
generally the case for the unpinned eigenvalues.) 

As before, in passing over to //; we restrict indices to 
enumerate different pairs. For N even, we find [consid- 
ering that the restricted sum has (N — 2)/2 terms] 



Sfii = \t - \tN\xi 



E 



Mi Mfc 



k^l,N+l 
= \t~ Tfl t - T(M? - 1 



Mi - Mfc 



E 



k^l,N+l 



1 

M; - Mfc : 



(A18) 



while if N is odd [where the restricted sum has (N— 3) /2 
terms] , 



S^i = - hrN^i - t 



k^tl,N,N+l 



MiMfc - 1 
Mi - Mft 



i T Mi 



r(Mf 



1) 



E 

kftl,N,N+l 



(A19) 



Mi - Mfc 



Furthermore, 



Ik- 



(A20) 



Comparison with Eqs. (All )-( |A13| ) shows that these 
are the same average increments, if we rescale r by a 
factor 2. The eigenvalues of UU* and O therefore exe- 
cute the same Brownian motion process, with the same 
stationary solution (A15). 



3. Brownian motion of orthogonal matrices 

As an independent demonstration of this correspon- 
dence between the eigenvalue distributions of UU* and 
O, we have investigated the Brownian motion of orthog- 
onal matrices. Let O be a random [N + 1) x (N + 1)- 
dimensional matrix in the orthogonal group, constrained 
to the sector Det O = — 1. 



4. Symmetry class BDI 

In class BDI we have X = U\WU*XU T , with U uni- 
form in U(N) and A a fixed diagonal matrix with entries 
±1 that sum up to Q. Since here the matrix X is sym- 
metric, X — X T , it is now diagonalized by an orthog- 
onal transformation; thus, the eigenvectors \k) = \k*) 
are real. As in class D, eigenvalues appear in complex- 
conjugate pairs, apart from eigenvalues pinned at 1. We 
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observe that fl mediates between the associated eigen- 
vector, \k) = S,k^\k) = ^ft*\k). In order to treat the 
partners symmetrically we have to require that that \k) 
is also real, so compensates any complex overall factor. 
It then follows that (k\£lfl*\k) = ^ = Xk, and thus the 
coefficients are related to the eigenvalues as in class 
D. 

To identify the pinned eigenvalues note that = W = 
f2 -1 is both Hermitian and unitary, and thus has eigen- 
values ±1. Let fl± be the eigenspace for each set of eigen- 
values, and fl± the analogous eigenspace for f2* , which 
is spanned by the complex-conjugated vectors. We de- 
note £ = signQ. The space [span(r2_^, fil^)]" 1 is then 
of dimension q = \Q\ (barring accidental degeneracies), 
and all of the vectors in this space obey X\k) = \k). 
Thus X has q = \Q\ eigenvalues pinned at 1. For 
each pinned eigenvalue, insisting that \k) = \k) implies 



il\k) — il*\k) — £|fc), £ — signQ = ±1 (consistent with 
the property that these states lie in the joint subspace of 
and 

With these additional properties in hand, the evalua- 
tion of drift and diffusion coefficients can proceed along 
the same steps as before. With the specified form of X, 
the incremental step of U carries over to an increment 

SX = iU[H, A][/ f rr - iQU*[H*,X]U T 

+ tQ{W + Q) + 2r(l - X) - 2NtX, (A21) 

where we already averaged terms of second order in H; in 
particular, terms such as UHXHWU* XU T = tQQ,* pro- 
duce the topological invariant Q. The associated eigen- 
value increment averages to 



5xi = - 2Ntxi + tQ{1\SI* + + 2r(l - x t ) 

- J^( x l ~ Xk)~ 1 (l\U[H,X]Wft* - HU*[H*,X}U T \k)(k\U[H,X}Wfl* - W*[H* , A] U T \ I) 

k 

= 2Nr Xl + 2rg%+ 2r(l - *,) - 4^' ^ - ^ - ^ ± ' T ' )/2 ± *'^ fc , (A22) 

, xi- x k 



where the 5ik term can be dropped because of the con- 
straint k I on the sum. Note how Q changes to 
q = \Q\ because of the sign of the matrix element in- 
volving pinned eigenvalues. 

Again we find that eigenvalues at unity remain pinned. 
For the other eigenvalues, we separate out from the sum 
the q eigenvalues that are pinned, and sum over the M — 
(N — q)/2 pairs of unpinned eigenvalues, 

Wi= - 2Nr Xl + 2r(l - a:,) - 2r 2 - ^ + Xj) 

xi - Xj 

A Xl A V^" X k + X k - 2x l 

-4rq --4txi]> , 

xi-1 xi + xj - x k - xj, 

(A23) 

where the double-prime again indicates the exclusion of 
the pinned eigenvalues. Furthermore, 



SxiSx m = 8r((5 te - S lm xf). (A24) 
For the quantities /i/ = (xi + xi)/2, this gives 

/x? - 1 



S t M = -2qr(m + 1) + 2r(l - 3 W ) - 4tJ^ 



SfXiSfim = 8r(l - fii)Si T 



- fik 



(A25) 
(A26) 



The stationarity condition (A14| is now fulfilled for 

M M 
k=l Km=l 

(A27) 

which gives the joint probability distribution of the 
twofold degenerate, unpinned eigenvalues (j, n (n = 
1,2,... M). 

5. Eigenvalue density 



The probability distributions (A15l and (A27) are 
both of the form 



/3 

m I : 



M M 

P(vi,V2, ■ ■ -Mm) oc H(l+Mfc) a ( 1- Mfe) 6 Yl Iw - ^ 

k=l l<m=l 

(A28) 

with (} = 2,a = 1/2, b=\Q\- 1/2 in class D and f3 = 1, 
a = 0, b = |Q|/2 - 1/2 in class BDI. These are called 
Jacobi distributions, because the eigenvalue density p(/i) 
can be written in terms of Jacobi polynomials |40) . 

For small N it is quicker to calculate the eigenvalues 
density by integrating out all /x„'s except a single one. 
Keep in mind that \Q\ of the /i n 's are pinned at unity, 
and that the N — \Q\ = 2M unpinned /x n 's are twofold 
degenerate. (The products in Eq. (A28) run only over 
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these M unpinned pairs.) The eigenvalue density 



(E"=i - Mn)) is then given by 
p(/i) = |Q|*(A«-l) + 2Mp(A*), 



p(p) = J d^i J dfi 2 

d/j M S(/j, - fii)P(pi, p 2 , ■ • ■ Mm)- 



(A29) 



The delta functions satisfy f_ 1 S(fx± l)dp = 1. The 
average conductance follows from the eigenvalue density 
according to Eq. 

(G) = T 2 G / dpp(/x)[l + (1 - Tf - 2(1 - Dm]" 1 . 

(A30) 

This gives the small- N results in Eq. ( 10 ) and Fig. [2j 
The large- TV limit ( 13 ) is obtained from an integral 



equation for the eigenvalue density in the Jacobi ensem- 
ble gang, 

M J ^ dfiptf) In | M - //| = -§(1 - 2/(3) lnp(p) 

- | ln(l + fi) - - ln(l - p) + C + G(l/M). 

(A31) 

The constant G is determined by the normalization 



dfipQi) = 1. 



(A32) 



The solution is 



Mp(p) 



M 



+ i(l - 2//3)[50u + 1) + 5(p - 1)] + 0(1/M), (A33) 



M = M + (a + b)/0 - |(1 - 2//3). 



(A34) 



Upon substitution of the values for a, 6, /3 in the two 
symmetry classes, and transforming back from p to p, we 



find 



1 



iV 



+ i(5(p - 1) - |<5(p + 1) + C(l/M), 



(A35) 

independent of Q and for both symmetry classes D and 
BDI. The corresponding result for the conductance is Eq. 
(13), to order 1/N if the limit N — > oo is taken at fixed 



6. Large-voltage limit 

For completeness we also give the derivation of the 
large-voltage limit (|9| of the average conductance. We 
need to evaluate 



(G) = T 2 G / <# p{<f>) 1 1 - (1 - T)e* | , (A36) 
Jo 

with p(</>) = (JJ n S(cj) — (j) n )) the density on the unit circle 
of the eigenvalues e 1 "^" of a unitary matrix Cl. 

In class D the matrix = U is uniformly distributed 
in U(N). This is the circular unitary ensemble (CUE, 
/3 = 2). In class BDI the chiral symmetry enforces that 
Cl is unitary symmetric, ft = UU T with U uniform in 
U{N). This is the circular orthogonal ensemble (COE, 
(3 = 1). Unlike the probability distributions we needed 
for the zero- voltage limit, these two distributions are in 
the literature 1401. 



JV 



■4>n)oc Yl 



J4>i 1/3 



(A37) 



k<l=l 



The corresponding density 



p(<P) = N/2ir, < < 2tt, 



(A38) 



is uniform irrespective of the value of (3 and without any 
finite- N corrections. Substitution into Eq. (A36) gives 
the result 
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